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ABSTRACT 

The accretion of hot, dense matter which consists of heavy nuclei, free nucleons, degener- 
ated electrons, photons and neutrinos is studied. The composition of free nucleons and their 
chemical potentials are provided through the equation of state for hot, dense matter proposed 
bv lLattimer & Swestvl ( I1991i ). The numbers of leptons are calculated through the reactions of 
neutrinos and through the simplified transfers of neutrinos in the accreting matter The ther- 
mally equilibrium fluid structure of the accretion is solved. The angular momentum transfer in 
the disk is analyzed with the shear stress which is driven in the frame of the general relativity. 
When the mass of a central black hole and the accretion rate are selected as Mbh ~ SM© 
and M w O.IMqScc^^, which provide the typical luminosity of gamma ray bursts(GRBs), 
the fraction of free protons in the accreting matter becomes very small, Yp « 10^**, while 
that of free neutrons is closer to unity, Yn ~ 0.7. Then the antielectron neutrinos v^. can 
freely escape through the disk but the electron neutrinos Vg are almost absorbed in the disk. 
Thus the frequent collisions of with Vg over the disk couldn't be occurred. The accretion 
disk is cooled mainly by Ve, which suppresses the increase of temperature and increases the 
density in the accreting matter such as T « 3 x 10^°K and p « 3 x lO^^g/cm'^ at the in- 
ner side of the disk. The scattering optical depths of and v^, then reach to be very large, 
Ts{ve) ~ Ts{ve) ~ 10^. Thus the accretion disk could be thermally unstable within the du- 
ration of diffusion time of neutrinos, idi// ~ lO(ms). The ram pressure produced by many 
scattering of Ve are very strong, which might produce the neutrino driven wind or jets. The 
luminosity and mean energy of neutrinos, L^, E^, ejected from the disk increases with the 
specific angular momentum of a black hole. The mean energy is inversely proportional to the 
central mass, Ei, cx M^^j, while the luminosity L^, is independent of its mass, Mbh- The 
diagram determining the physical properties of a central black hole from the observation of 
each flavor of neutrino is proposed. 

Key words: accretion, accretion disks — black hole physics — gamma rays: bursts — neu- 
trinos. 



1 INTRODUCTION 

An intense gamma-ray burst, GRB030329, allowed us to carry 
out the very detailed observations of optical, x-ray, and radio 
counterparts to the burst. Its photometry and spectroscopy led 
to the separation between the afterglow and SN contribution s 
jMatheson et al."2003': Kawabata et al. '2003"; 'Lipki n et alj|2004) . 
The SN component is similar to SN 1998bw; an unusual Type Ic 
SN fstathakis et alj|2000l: iFerdinado et alj|200lh . It makes clear 
that at least some GRBs arise from core-collapse SNe. GRB030329 
was seen as a double-pulsed burst with the total duration time 
of ~30 sec JVanderspek et alj |2004'). Its total isotropic energy is 
E~i ~ 10^^ erg. The GRBs with known redshifts distribute the 
emission energy over the wide range of energy such as 10** 
erg ^ ^ 10^'*erg jGreine3l2004] . The duration times of GRBs 
extend from ~IOms to ~ 10^ sec lPiran 1999 : Stem et alJl200ll: 
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iLlovd-Ronning & Ramirez-Ruizll20o3) . The light curves of GRBs 
show really various profiles. However, to fully understand the emis- 
sion mechanisms of the observed 7-ray bursts and of the afterglows 
one must have a firm handle on the generation process of the 7-ray 
bust. 

The fireball model of the 7-ray bust has succeeded in un- 
derstanding the characteristics of some light profiles of GRBs 
jPirarJ Il999l) . Its model was proposed to solve the compact- 
ness problem, that is, to understand both the short time vari- 
ation of 7-ray radiation and the optical ly thin process in the 
emitt i ng region of the 7-ra v photons dGuibert. Fabian. & Reed 
Il983t ICarigan. & Kat3ll992h . The solution of the fireball model 
requires the beaming jets with the emitting region relativis- 
tically moving to our direction. The two processes are pro- 
posed for the formation of the jets, one is the magneto- 
hydrodynamical process Thom pson. Chang. & Ouataerll l2004t 
ILvutikov & Blandfor j2003 ), and the other is the neutrino process 
TPopham. Wooslev & Frveill999tllNaravan. Piran. & Kumatl200lt 
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iKohri. & Mineshigd l2002t iDi Matteo. Pema. & NaravanI l2002h . 
Further more the two different roles are expected for the neutrinos 
to form the jets. The frequent colhsions between neutrinos and an- 
tineutrinos produce the fireball with light elements which expands 
relativistically in a funnel around a rotation axis. The other role is 
the acceleration of plasmas by the collisions of neutrinos with nu- 
cleons or electrons, which produces the neutrino driven jets. The 
sufficient collisions between u and u at the central region of the ac- 
cretion disk were proposed by .Popham. Woosley & Fryer ( 1999i) . 
They assumed t he accretion disk is optically thin for the transfer 
of neutrinos. But. lPi Matteo. Perna. & NaravanI IzOOA showed that 
the accretion disk becomes optically thick for neutrino transfer and 
then the neutrinos trapped in the matter fall into a black hole. 

However, the composition of accreting matter should be 
carefully considered. The production rates of neutrinos and 
their opacities directly depend on the density of free nucleons. 
[Kohri, & Mineshige (2002) were thought in question whether the 
neutrinos can cool the accretion disk. They treated the accreting 
material in /3 equilibrium and showed that there exist just a small 
number of free nucleons and then the neutrino processes are no 
longer effective in cooling of the disk. This result is in a striking 
contrast to the previous works on the accretion disk in which the 
material almost consists of free nucleons produced by the photo- 
disintegration process and the dominant cooling was carried out by 
neutrinos. Material at high temperature T ^ 5 x lO^K becomes 
in nuclear statistical equilibrium. Its composition of material de- 
pends on the number ratio of photon to baryon, i/; = / hhtib, 
where jiH is the mean molecular weig ht, n-y and tlb are the num- 
ber densities of photons and baryons ^ Meveretal]|l992h . If ■(/> is 
large, tp ^ 1, the baryons form a dissociated gas of nucleons 
and alpha particles. However, if ^ is of order unity or less, iron- 
group nuclei are typically favored. The ratio tjj is expressed as 
■tp — 0.034^^1 /pi2 in terms of the temperature Tn in units of 
lO^^K and the density pi2 in units of lO^^g/cm'^. In the case of a 
typical accretion with M = O.lM0/sec and Mbh = SMq, the 
temperature and the density become are 10^" < T < 10^^ (K) and 
10^^ < p < 10^'^(g/cm'^) at the inner region of the accretion disk. 
Then the values of V' become as i/) = 3.4 x 10~^~~"^. One should 
treat the accreting matter with a large abundance of heavy nuclei. 
We determine the densities of free neutron and proton, alpha parti- 
cle, and heavy nuclei by using the equatio n of state for hot, dense 
matter given bv lLattimer & Swestvl il99 it) . The changes of lepton 
number per a baryon, Vj, are calculated through the reactions of 
neutrinos with these baryons. 

The thermal equilibrium structure of the accretion disk 
should be determined in th e frame of a general relativity. 
iMacFadven & WooslevI Jl999h showed that the core-collapse in 
the evolved star with the main sequence mass , M 35M0, 
forms a black hole at the center with the rotating disk balanced 
with gravity when the evolved star has the some large specific an- 
gular momentum. Then the specific angular momentum parame- 
ter a of the falling matter into a black hole becomes larger to 
be a ^ 0.9. The falling velocity of the rotating disk is re- 
quired to be very small in comparison with light velocity, since 
the time duration of some GRBs is a few seconds, i.e., viaii ~ 
scaleofdisk/durationtime « W^rg/afewsec ^ 10~^c, 
where Vg is the gravitational radius of a black hole and c is the 
light velocity. The accretion disk with slow falling velocity should 
be investigated in the frame of rotating spacetime. The heating rate 
due to viscosity is proportional to the shear stress anp. In the static 
spacetime its stress increases with the differential rotation of the 
accreting flow a^p oc drO,, where SI is the angular velocity of the 



falling flow. In the rotating spacetime the differential rotation of 
the fluid observed in the frame co-rotating with spacet ime is phys- 
ically meaning ful I'Thorn e. Price. & Macdoniij3ll98d) . In the ex- 
treme Kerr the relative angular velocity Q — uj has the maximum at 
the radius r ~ 2rg when Q is the angular velocity of the Keplerian 
orbit and uj is the angular velocity of spacetime. It is convenient 
for the analysis of energy balance and of the transfer of angular 
momentum to introduce the orbiting frame with orthogonal tetrad 
which rotates with Keplarian velocity. In this paper the shear stress 
ar^p and the transfer of angular momentum are solved by using the 
orbiting frame. 

The power by advection contributes seriously to the energy 
balance in the accretion disk. The previous works on the accre- 
tion disk with neutrino loss had be en studied in the frame of an 
advection-dominated accretion fl ow tDLMatteOjJema,_&>iaraYan| 
200?; 'Kohri. & Mineshiae '2002!; iNaravan. Piran. & KumaJl200lt 
Ponham. Woosley & Fryer 1999). The role of the advection in the 
energy balance depends on the velocity profile since the power by 
the advection is expressed as PV-v, where P and v are the pressure 
and velocity of the flow. When V ■ w > 0, the accreting matter is 
cooled by the expansion of fluid. On the other hand when V ■ « < 0, 
it is heated by the compression like as in the core-collapse in a su- 
pernova explosion. The radial velocity profile of the accreting flow 
is determined by the energy balance of the disk and by the transfer 
of angular momentum. The scale height at the inner portion of the 
disk with the temperature T = lO'^'^^^^K is very thin. The angular 
velocity of the accretion disk is then approximated to be Keplerian. 
We determine the radial profile of the velocity by solving the trans- 
fer of angular momentum with the thermal energy balance of the 
disk and evaluate the efficiency of advection in the energy balance. 

The thermal stability of the accretion disk might depend on 
the scattering optical depth of i/e- The thermal energy in the disk 
is carried out mainly by Pg. The collision frequency of i>e with the 
nucleons or electrons reaches to N ~ r|^^ ^ 10* within the disk 
for the typical case of the accretion. The escaping time from the 
disk is then tesc ~ Ta^aH/c ~ lO(ms), where H is the height 
of the disk. The thermal energy produced by viscous heating stays 
within the disk in the time At « tesc and might induce the thermal 
instability. We shall briefly analyze the thermal instability. 

If the luminosity of each flavor of neutrino and its mean energy 
are observed the physics of GRBs will be cleared. We shall calcu- 
late the luminosity L^. and the mean energy E,^. for each type of 
neutrino i/i emitted from the accretion disk. The diagram determin- 
ing the mass of a central black hole, Mbh, and its specific angular 
momentum, a, from the observations of neutrinos is presented. 

In this paper we study the accretion of hot dense matter onto 
a black hole. In Sect. 2 we describe the equation of state for the 
accreting hot, dense matter. The reactions of neutrinos and the 
changes of neutrino fractions are presented in section 2.2. The rel- 
ativistic model of accretion disk with neutrino loss is given in Sect. 
3. In Sect. 4 the results of thermally equilibrium disk are given. The 
profiles of compositions along the accretion flow are shown in sec- 
tion 4. 1 . The flow structure and chemical potentials of nucleons are 
depicted in section 4.2. The emissivities and opacities of neutrinos 
are presented in section 4.3. The luminosity and mean energy of 
neutrinos versus mass and angular momentum of a black hole are 
discussed in section 4.4. In Sect. 5 the thermal stability of the accre- 
tion disk is investigate. The neutrino transfer in homogeneous disk, 
changing rates of positrons and the dynamical properties of accre- 
tion disk are presented in Appendix. The cross sections of neutri- 
nos, opacities, emission rates and the reaction rates of neutrinos are 
summarized in Tables, which are shown in Appendix. 



Hot, Dense Matter Accretion onto a Black Hole 3 



2 ACCRETING DENSE MATTER AND REACTION OF 
NEUTRINOS 

2.1 Equation of State for Accreting Dense Matter 

It is important for the study of thermal and dynamical properties of 
massive accretion disks precisely to determine the fraction of free 
nucleons in the accreting matter Its fraction provides the emissivi- 
ties and absorptions of neutrinos. The mass fraction of free n ucle- 
ons depends on th e photon-to-baryon ratio (A cBurbidee et aljl957t 
lMeveretalJl992h . 

(1) 



^ pNa' 

where Na is Avogadro's number We compute in terms of the 
temperature T and density p as 



37 C(3)(fcr) 



7r2 (/ic/27r)3 pNa 



3 rpZ 
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where g-y is the helicity for photons. We evaluate the ratio (j> 
for the neutrino dominated accre ting flow (NDAF) given by 
iDi Matteo. Pema. & NaravanI <2002h . The loci of NDAF with M = 
O.lM0sec~^ and the contours of (f>{p,T) in the p — T plane are 
shown in Fig 1. The ratio (j) on the flow is </> 0.1, showing 
little number density of photons. On the othe r hand the previous 
works on NDAF without iKohri. & Mineshigd llOO'j) adopted the 
approximate formula for the mass fraction of free nucleons, Xnuc, 
(Oian& Wooslev 1996) : 



9/8 



■exp(-0.61/rii) 



(3) 



The contour of Xnuc is also plotted in Fig. 1. At the inner side 
of NDAF its fraction, Xnuc, reaches to unity. We also calculated 
the flow profile of the standard accretion disk with the Keplerian 
angular momentum in which the fraction of free nucleons is pro- 
vided through this formula, Xnuc- The accreting velocity of the 
standard disk is much less than that of NDAF, which produces the 
denser profile of the accretion flow for the same accretion rate, 
M — O.IMqscc"^. Then its ratio cj) is very minute, 4> ^ 0.01. 
Thus it is necessary to handle the accreting matter with an appre- 
ciable abundance of heavy nuclei. 

We adopt a simplified, analytical equation of state for hot, 
dense matter given by .Lattimer & Swestv ( 1991.') . in which the 
heavy nuclei are treated as being in a b.c.c. lattice. In accordance 
with the Wigner-Seitz approximation each heavy ion is consid- 
ered to be surrounded by a charge-neutral spherical cell consisting 
of a less dense vapor of neutrons, protons and alpha particles as 
well as electrons. The volume of this cell is given by 14 = , 
where ua is the number density of heavy nuclei. We construct 
a three-dimensional table of the relevant thermodynamic quanti- 
ties, i.e., nn,np,na, p.n — Pp, A, Z, Pt, as a function of the inputs 
(p, T, Yc), where the number densities of free neu- 

tron, free proton, a particle, pn — Mp is the difference of chemical 
potential between the neutron and the proton, A and Z are the mean 
mass number and proton number of heavy nuclei. The contours of 
chemical potential, pn — Pp, are depicted in Fig.9. 



2.2 Reactions of neutrinos 

The thermal and dynamical properties of a massive accretion disk 
are greatly influenced by neutrinos. The neutrino loss cools the 
disk and the interaction with ambient matter may produce the con- 
vection. In this paper the dominant neutrino reactions in the stel- 
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Figure 1. The loci of accreting gas on the p(g cm"'') — T(K) plane. NDAF 
is the locus of the neutrino-dominant accretion flow with M = O.IMq 
sec~^ and Mgu = 3Mq. SADl is the locus of the standard accretion 
disk with Keplarian angular momentum where heavy nuclei are ignored. 
SAD2 is the locus of the standard accretion disk with heavy nuclei. The 
lines with the constant photon-to-baryon ratio ip'yb and those with constant 
mass fraction of free nucleons Xnuc are plotted. The phase boundary be- 
tween heavy nuclei and uniform evaporated matter for the case Yg^ = 0.3 is 
described by the dot-dashed line with the label Pnp,A = Pnp,G- 



lar core collapse are included, which are listed in Table 1. The 
summaries of the various ne utrino cross sections of relevance in 
supern ova theory are given in lTubbs. & Schramml il975h . lBruennl 
^198^ and lBurrowd j200ll) . The approximate blockings in the lep- 
ton phase space are described in Ruffert, Janka, & Schafei ( 199q). 
Here the leptons spectra are assumed to be represented by ther- 
mally equilibrium Fermi-Dirac distributions with the temperature, 
T (taken equal to the gas temperature), and the degeneracy fac- 
tor, rji. The neutrino cross sections, opacities, emissivities and re- 
action rates used in our model are listed in Table Al-5(Appendix 
A). When the density ni and temperature T of leptons are given, 
their chemical potential pi is determined through the Fermi inte- 
gral, ni{rii) = AnikT/hcf dxx^/[l + exp{x - -qi)]. The 
electron-positron pair is assumed to be in equilibrium with thermal 
photons. 

The reaction and transfer of neutrinos are evaluated in the 
simplified disk with the thickness H{r) in which the matter is as- 
sumed to be homogeneously distributed in the vertical direction. 
The thickness H{r) is taken approximately to be a scale height. 
The change of lepton fraction, Y, is calculated along the accret- 
ing flow. At the inner region of the disk, the hot, dense matter 
changes the lepton fraction within the duration time of the matter 
moving along the distance Ar with the accreting velocity Vaccrt, 
At = Ar/vaccrt. When Vaccrt ~ 10~''c ~ 10^ [cm/sec] and the 
section of the moving distance is taken to be Ar ~ O.lrms ~ 
lO''(Msjf/3M0)[cm], the duration time is At ^ 10[ms]. In this 
time At the neutrino reactions do not always reach to the local 
thermodynamic equilibrium(LTE). The changes of leptons are cul- 
culated without LTE. 

The changing rate of each lepton number in unit volume, ni , 
is produced through the networks of reactions, 

he = -hp + h^-riA + h'^iiCA, (4) 
= —hjj ~ hce, (5) 

hi,^ = -I- + nee + njvjv + njiicyi + 'i-7 , (6) 
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Table 1. Netrino reactions 



No. 


Neuti'ino reaction 


Reaction process 


References 


1 


e4 


- p i/e + n 


/3 process 


Burrows(2001) 


2 


e+ 


+ n ^ Ue + p 


positron capture on neutron 


Burrows(2001) 


3 


e4 




electron capture on nuclei 


Bi-uenn(1985) 


4 


- 


-*v, + N 


neutrino-nucleon scattering 


Burrows(2001) 


5 


i^i - 


^A^v^ + A 


neutrino-nucleus scattering 


Bruenn(1985) 


6 


t^i - 


-\- e ^ Ui + e 


neutrino electron scattering 


Burrows(2001) 


7 


e4 


- e+ ^ Vi + Vi 


electron pair process 


Burrows(2001) 


8 


N ■ 


JrN^N + N + Ui + v, 


nucleon-nucleon bremstrahlung 


Raffelt(2001) 


9 


riH 


-n^n + p + e + Ue 


Urea process 


Shapiro and Teukolsky(1983) 


10 


7 - 


-+ j/j + pj 


plasma decay 


Ruffertetal.(1996) 



"-i^e = n-^ + riee + riNN + 
hij^ = hee + riNN + nl^ ■ 



(V) 
(8) 



Here n0,hfj,hEe,nA and n jv jv are the net production rates of neu- 
trinos by the following reactions; the /? process, e + p ^ Ue + n, 
the positron capture on neutron, e+ + n ^ + p, the electron 
pair annihilation, e + e'^ ^ Ui + Vi, the electron capture on the 
nuclei, e + A ^ + A' and the nucleon-nucleon bremstrahlung, 
N + N ^ N -\- N -\- Ui -\- Vi. The positive sign in these net rates 
indicates the increase of neutrinos. The plus sign at the upper suffix 
denotes the production rates of neutrinos without inverse reaction, 
i.e., n^nQj^ and are the production rates by the URCA process, 
N + N —> n + p + e + Ve, and by the plasma decay, 7 ^ /^i + Ui. 

The changes of neutrino density are also brought about by the 
leakage of neutrinos from the disk. The flux number density of each 
neutrino at the surface of the disk, F^j^,.), are approximately ex- 
pressed by the production rate of neutrino, fi^^ , and the absorptive, 
scattering and total opacities, k^.s and n^. = n^.a + Huis, 



n+H 



(9) 



where and r^.a are the total and absorptive optical depths for i^i 
(see Appendix B). The absorptive and scattering opacities for each 
type of neutrino are given by 

a s . s . 3 r 3 1 1 W 

l^u^a — N : ^^i3 — ' i ~r : t t J 

where Kfi,Kp, k\ and are the scattering opacities due to the scat- 
tering with neutron, proton, nucleus and electron which are repre- 
sented in Table A2. 

Thus the change of the each neutrino fraction, . , along the 
accreting flow is expressed by 

ribYv- = II Vrn^. = ^7^. (12) 



H • 

The changes of the electron fraction. Ye, and of the positron frac- 
tion, Ye+ , are given by 



nbYe 



K+ = 



ne+th 



(13) 



Uf, Ub Ub 

where n^+^f^ is the number density of positrons which is equilib- 
rium with thermal photons(see Appendix C). 

The energy loss carried out by neutrinos is similarly expressed 
with the emissivities of neutrinos, . . The emissivity of each type 
of neutrino is described as 



Qv^ = qi3 + qee + QA + QNN + qu RC A + q-, , 

<li^e = 10 + Qee + gjVJV + qj , 



(14) 
(15) 



qeeil^x) + qMN^Vx) + q-fi^x), 



(16) 



where qp,q0,qee, qA,qNN , quRCA and q^ are emissivities by the 
f3 process, the positron capture on neutron, the electron pair anni- 
hiration, the electron capture on the nuclei, the nucleon-nucleon 
bremstrahlung, the URCA process, and by the plasma decay; 
qee{i'x),qNN{i^x),q-y{i'x) are those for heavy neutrinos, Ux = 
Uf_i,Vr (see Table A3). The energy flux density of each neutrino 
at the surface of the disk, , is then given by 



(17) 



3 RELATIVISTIC MODEL OF ACCRETION DISK WITH 
NEUTRINO LOSS 

The duration time of GRBs requires that the collapsed, accreting 
matter has large specific angular momentum. The central black 
hole may be rapidly rotatingi MacFad ven & Woosl ev 1999). We 
derive precisely the shear stress tensor cr^r in Kerr spacetime. 
The relativistic model for the accretion disk was investigated by 
Novikov, & Thome 1 1973). However, their model is insufficient in 
some respects. The boundary condition at the inner disk is not suit- 
able, i.e., the thickness of the disk becomes infinitesimal and the 
density diverges to infinity at the boundary. The work produced 
by advection was not included in their model. We reconstruct the 
relativistic model without singularity and with advectional motion 
based on the standard theory of the thin accretion disk. 

We choose units with G = c = 1. The metric of spacetime, 
in terms of Boyer-Lindquist time t and any arbitrary spatial coordi- 
nates X-' , has the form 



ds^ = -adr + gjk{dx^ + P^dt){dx'' + P^dt). 



(18) 



The metric coefficients a, {3^ , gjk are the lapse, shift, and 3-metric 
functions. Theses functions are given by 



[ A ) ' 



9r^ ^ ' 



^1 



sin^ OA 



13" 



-u — — ^"^^""^ ^ pi — g.^ — Q other j and k, 



where M is the mass of the black hole, a is its angular momentum 
per unit mass, and the functions A, E, A are defined by 



A = r^-2Mr + a^, E = 4- cos^ 6*, 



(19) 



We introduce a set of local observers who rotate with a Kep- 
lerian, circular orbit. Each observer carries an orthonormal tetrad. 
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For the Keplerian orbiting observer with a coordinate angular ve- 
locity of Jl, its world line is r ^constant, 9 =constant, = 
constant. The observer in the LNRF who rotates with the angular 
velocity uj measures the linear velocity of a particle moving in a 
Keplerian orbit as W(^) = ^g^^{Q, — lo). Its Lorentz factor is 



7 = (1- 



-1/2 



We adopt the Keplerian orbiting, observer's 



frame ("orbiting frame") with the set of it's basis vectors: 
eg = 76' 



/ 9 ^ d \ 



e,-. — e 
where e'\ , 



d_ 

dr ' 



dd' 



(20) 
(21) 
(22) 



2M2 



are difined by e" = a, e^"^ = g^^, e^''^ = 



The corresponding orthogonal basis of one-forms is 

!* = adt - •!;(^)7e*dv?, Cj'^ = -^-ye^At + -ye^dtp 



(23) 



where a is the "lapse function" in an orbital frame, a — 07(1 -I- 
The shear stress in an orbital frame is expressed by 



(24) 



where a°'^ is the shear stresses in the Boyer-Lindquist coordinate, 
and is the Boyer-Lindquist components of the basis of one- 
forms 1123 1 in orbital frame, cZ5" = e'^Ax°' . The shear stress in the 
Boyer-Lindquist coordinate is defined by 



(25) 



where /i^^ is the projection tensor, h^v = Qtm + u^Uy. The shear 
stress in the orbiting frame is then represented as 

d\n(^ — uj) 1 dlnu) ^ 



2 L a/ Qrr 



dr 



1 



dr 



(26) 



where is the rotating velocity of spacetime, — a~ ^g^^uj. 
While the absolute value of the shear stress in the Shwarzchield 
spacetime monotonously increases as |a'^'^| « 0.7511 according to 
r rms, the shear stress in rapidly rotating spacetime has the 
minimum in the vicinity of rms- 

The generation rate of the viscous heating is evaluated in the 
orbiting frame by 



(27) 



where f^- is the stress tensor. For the "a model" of the viscosity 
the stress tensor is expressed by t-.^ — a^iaPth, where Pth is the 
thermal pressure and a-uis is a parameter with a^is = 0.1 — 0.01. 

Then we consider the energy equation in the orbiting frame. 
The energy-momentum t ensor of the viscous fluid with heat flux 
is express edl Misner, Tho ne & Wheelei)ll973t iNovikov. & Thorn3 



IS exp i 
Il973t f 



lYokosawa.1995.) as follows, 



S/3 



= (p + e)u u + Ph ~\- 1 +g u -\- u q 
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(28) 



where p, e, and P are the rest-mass density, internal energy den- 
sity and pressure and q" is the heat flux carried by neutrinos and 
photons. By using the covariant derivative, V, the local energy con- 
servation, it • (V ■ T) = 0, is represented as 



+ _V; ■ a{eu + q') + P r 



(29) 



The mass conservation V ■ (pu) = is written by 



(30) 



+ ^V.apn^ = 0. (31) 

By using the projection vector in ip direction, h''' = h''''^, the Euler 
equation in the ip direction, h''' • (V • T) = 0, is expressed in the 
frame fixed at the distant stars. 



di 



+ -\/i(a(p + e + P)u*u^+tt+u^q*) (32) 

dap 



+ - 



idxf 



= 0. (33) 



We are interesting in the gross properties of the disk with neu- 
trino loss. For simplicity, we consider a steady state disk with axi- 
ally symmetry. As in the standard theory of thin accretion disks, we 
consider height-averaged quantities. Putting the integrated stress as 
W = tftpdz — OvisPthdz, and introducing the flux den- 
sity at the surface of the disk, F = q'^ (z = H), the energy equation 
is expressed by 

F + (a,:^ + {\/r&u^ + P-Wfaeu^)) W = 0. (34) 

\ aavis / 

From the equation <33l l the conservation of angular momen- 
tum is 



Vi ■ (a(p + e + P)uu^ + + u^q^ — 



(35) 



The angular momentum is transported by the viscous stress and 
by the neutrino flowing g'. Introducing the total mass flux A/, we 
express the transport equation of angular momentum as 



dr^ 



M 
— 7 

2ti 



It has the solution: 
M, 



W - 



2-K 



nil, 



(36) 



(37) 



where II is the function with the value of unity at the great distance 
from the hole: 



n 



u^^C ^ 2-Kfg 



(38) 



The term }g is due to the neutrino flowing. The functions / and ; 
are 



f = 



M 

'2^ 



(Tfip^l + ( 



2vju^Of^ 
exp 



{u^ + C)dzu, 



2zuUu,af0 

■dzu 



1 + 



P 2(jf.0avisTAj- 



dvj 



+ 



aiiK 



(39) 

(40) 

(41) 
(42) 



l + P/e dzu 

Here C is the integral constant. Novikov & Thome(1973) adopted 
such as C = —Uip{rms) since they assumed no viscous stress can 
act at r = rms - The accreting matter with viscosity co ntinuously 
flows into a black hole through the boundary radius r.r.. iYokosawal 
I l99^ . We set the constant C so that the accreting matter continu- 
ously flows through the radius rms and then the viscous stress has 
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the value such as tr^{rms) = Qi,isp(rms), 
47r/i 



C = -Mv(rms) + ^{rms)^(rms)e* {rms)0lvrsP{rms)- (43) 

Dynamically equilibrium structure of the disk is given by the 
Euler equation. The equation of motion of the fluid, h^T^.j = 0, 
for the steady, axially symmetric and rotating fluid can be expressed 
in the total differential form: 



-dU. 



(44) 



p + P 

The equipressure surfaces are given by the equation U =constant. 
The quantity U = U{p) is equal in the Newtonian limit to the total 
pot entieJ (gravitational plus centrifu gal) expressed in the units of 
^Abra mowicz. Jaroszvnski Sikora 1978). In the neighborhood 
of equatorial plane the potential U is expanded with z: 



U = Uo{r,0) + -nlz\ 



where Qq is 



GM 



$(r) 



(45) 



(46) 



The function $(r) becomes unity far a way from a hole. Its explicit 
expression of ^(r) is shown in Appendix D. The thickness of the 
disk is approximately given by the scale height, 



H 



(47) 



The energy equation 1341 integrated over the thickness of 
the disk is represented with the mass accretion rate, M = 



F + (Jr^W + 



M 



d e + P PdlnP 



= 



(48) 



p p dlnr , 

The flux density F is evaluated by the simplified neutrino 
transfer, 



F 



(49) 



Ve, V^, V^i, ffii ^T, > 



where Vi means the flaver of neutrino, Vi 
Thus the basic equations J48t and 13711 determine the stationary, 
thermally equilibrium structure of accretion disk. 

In addition to the above formalism of accretion disk, we shall 
consider the thermal energy of the partially degenerated matter 
which contributes to the viscous heating and the advection. In the 
ranges of density and temperature, 10^^ < p < 10^''(g/cm'') and 
10^" < T < 10^^ (K), the degeneracy factors of electrons has the 
values in the range, 5 < 77,; < 15 (see Fig. 10). The pressure and 
the internal energy of electrons are expressed as 



P. = 



Me 

127r2 
kT 



1 



7 2 
— n 
15 



Pe 
47J-2 



fcT 

Pe 
1 

" 2 
kT 

Pe 



27r" 



kT 



kT 



(50) 



kT 



+ 



kT 

Pe 



7 2 
15" 



kT 



(51) 



The first terms in the both square brackets provide the pressure 



and internal energy corresponding to the completely degenerated 
state. The second and third terms depend on the temperature. We 
introduce the thermal pressure of electrons, Peth, expressed by the 
above expression without the first term. The completely degener- 
ate matter has no contribution in cooling or heating by advectional 
motion. Thus the pressure and internal energy density without the 
completely degenerate parts are used for the advection or for the a 
model of viscosity. 

The baryonic matter is hardly degenerated over the above 
phase area in p — T plane. The pressure of baryonic matter is con- 
sist of that of free nucleons P„uc, that of a particles Pa, that of 
heavy nuclei (Lattimer & Swestv 1991), 

Pb = P^uc + Pa~ P(D~uD') + ^h[kT{l-u)~ pAu],{52) 

where pn is the chemical potential of heavy nuclei, u is the fraction 
of space occupied by heavy nuclei, f), D, D' , h are the functions of 
nuclear quantities iLattimer & SwestMll99lh . The values of Ph in 
the above phase area are almost same as ni,kT. Thus the internal 
energy of baryonic matter is approximately set to be et ~ 3/2Pf,. 
The total pressure and internal energy density are given by 



P = 



4,11 21 
e = eb + €e + aT [— + —t), 



(53) 
(54) 



where a is the radiation constant r is the neutrino opacity, r — 
E(r^j2 + 1/V3)/{t^J2 + l/VS + l/3r^,a). The third terms 
on the right hand side represent the contributions of radiation in 
thermally equilibrium with relativistic electron positron pairs and 
of neutrinos lPopham. & Naravan 1995.) 



4 RESULTS 

We have solved the stationary flow structure at the position range, 
r — rm.s ~ 10^r„is. At the outer boundary of the calculation we 
select the most large value of Ye within Ye ^5 0.5 which satisfies the 
stationary conditions, <48> and <37t . The fractions of neutrinos at 
the boundary, Y^^ , are taken which provide the local thermodynam- 
ical equilibrium in neutrino reactions. The accretion rate M is taken 
to be Af — 0.01 ~ IMqscc"'^, and the angular momentum param- 
eter of a black hole ranges from a = ~ 1. The accreting velocity, 
V = w'^/u", is derived from the mass conservation, AI — const., 
for given density p and height H. The typical parameters in the 
models of GRBs are selected as a = 0.9, M — O.lMosec^^, and 
Mbh = SMQjMacFadven & Wooslev 199^. The viscous param- 
eter is set to be ct^is = 0.05 jTumej2004l) . 



4.1 Flow structures 

The flow profiles are shown in Fig. 2, where the density p, temper- 
ature T, accreting velocity v and scale height H are plotted in the 
radius normalized by the gravitational radius, r , = r/rg. While the 
density p increases according to the accretion rate M and reaches 
to a very high value, p ^ 10^^[g/cm^] for M ^ O.IM© sec , the 
temperature, scale height and accreting velocity are little changed 
for the wide range of accretion rate, M = 10~^ ~ IMqscc^^. 
The profiles of temperature and accreting velocity are nearly con- 
stant along the flow, T 3 x 10"'K and v ^ 10"^c. The scale 
height is expressed as « 0.1rg(r,)^'^*, showing the very thin 
disk formed around a black hole. 
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Figure 2. The flow profiles of density p(g cm~"^), temperature T(K), 
drift velocity i)(cm/sec), scale hight //(cm). The radius r is normalized 
by the gravitational radius Vg. The mass of a black hole is fixed to be 
Mbh = 3M0, and the angular momentum parameters are selected 
such as a = 0, 0.9, 1. The numbers denote the accretion rates, M= 
1,O.5,O.2,O.1,O.O1(M0 sec-l). The flow profiles with M = 0.1M©sec-l 
are plotted by filled lines, and those with M = IMqscc"^ and M = 
lO^'^MQsec"^ are depicted by dot-dashed lines. 
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Figure 3. The same as Fig. 2, but for the mass dependence of a black hole. 
The numbers denote the masses of a black hole, Mgu= 1,2,3,5,1O(M0). 
The accretion rate and the specific anguler momentum of a black hole are 
fixed as M = O.IMqscc^^ and a=0. The radius r is normalized by the 
gravitational radius rg . 

These flow profiles are somewhat different from the pre- 
vious works on the advection dominant flow with neutrino 
lossCP opham. Wooslev & Fryer 1999; Narav an. Piran. & Kumar 
l2001t iDi Matteo. Perna. & NaravanI 12002) . In the previous cases 
the temperature increases in proportion to the accretion rate. The 
profiles of the temperature and velocity rapidly increase to be 
T ^ 10 K and « — > c for r — > Vms- Then the density does 
not reach to the extremly high value which is restricted to be 
p « lO^'^^^'^ [g/cm'^] in the ranges of accretion rate, m = 0.1 ~ 
lOM0/sec. These differences in our results and previous ones are 
caused mainly by the cooling process, which will be discussed later. 
In the previous cases the most of neutrinos are absorbed by free nu- 
cleons and thus trapped in the accreting flow. In our study the anti- 
electron neutrinos Ve are little absorbed by free protons since free 
protons are reduced in the fraction, Yp ^ 10~^. The accretion disk 
is then efficiently cooled by Ue- It suppresses the rise of temperature 
and velocity, which then produces the highly dense flows. 



-1 




-7 ' ' ' • • ' ' ' 
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Figure 4. The profiles of total pressure P, electron pressure Pe, pressure 
of baryonic matter P;,, radiation pressure Pr, total thermal pressure Pth, 
and the thermal pressure of electrons P^th are plotted in the case of the 
accretion with M = O.IMq sec'^, Mbh = 3Mq and a = 0.9. The 
unit of pressure P is (MeVfm-3). 

The density p and the temperatute T decrease as the mass of 
a central black hole increases when the accretion rate is restricted 
to be constant. The flow profiles are shown in Fig. 3 where the 
mass Mbh is changed from IMq to IOMq and other parame- 
ters are fixed to be a = 0, A4 = O.lM0sec~^. Their scale de- 

r /3 

pendences are described approximately as T oc Mgj^ , H oc 
Mbh and p oc MMg'jj. That of the surface density is then 
E = pH oc MgljM. Thus the obtained flow profiles are ex- 
pressed approximately by p « 2 x 10^*r-~^'^mm~^[g/cm^] and 
T ~ 5 X W^^rZ^'''^m~^^'^K, where m is the mass of a central 
black hole normalized by IMq and m is the accretion rate normal- 
ized by lM0/sec. 

The inner disk is supported mainly by the baryon pressure Pt- 
The compositions of the pressure, Pb,Pc, Pr and the thermal pres- 
sure, Pth, Pe,th are shown in Fig. 4. While the electron pressure 
Pe dominates at the outer region, r ^ ISr^s, it more gradually 
increases than Pb according to the accreting flow. When Ye rapidly 
decreases, the electron pressure itself decreases along the flow. The 
contribution of the radiation pressure to the total pressure is negli- 
gibly small at the inner region, but at the outer region, r ^ lO'^rg, it 
becomes comparable with that of Pb. The thermal pressure of elec- 
trons Peth is comparable with its degenerate part of the pressure. 
The degeneracy factor of electron rje is ranged as 77e ~ 2 ~ 10 for 
m = 0.1 which doesn't reach to an extremely large value(see Fig. 
10). 

4.2 Composition and chemical potential 

The composition of accreting matter which consists of free neu- 
trino, free proton, electron, positron, and each flavor type of neu- 
trino is shown in Fig. 5-8. The massive accretion rates, M — 
0.05 ~ lM0sec~^ for a — 0, are selected in Fig. 5 and Fig. 6, 
while the less ones, M = 0.005 ~ O.IMqscc"^ for a = 0.9, 
are shown in Fig. 7 and Fig. 8. The fraction of free neutron Y„ in- 
creases along the accreting flow and reaches to a dominant value 
at the inner side of the disk, F„ ~ 0.7 ~ 0.8. However the frac- 
tion of free proton decreases along the flow and reduces to an ex- 
tremely small value, Yp ~ 10^'^^^*. When the accretion rate m 
increases, the proton fraction decreases as, Yp oc though the 
fraction of neutron remains a constant value at the inner side of 
the disk, Yn ~ 0.7 ~ 0.8. The number ratio of free neutron 
to free proton, iin/np, is proportional to the degeneracy factor. 
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Figure 5. The profiles of electron fraction Ye, positron fraction Yj,+ , free 
neutron fraction Yn, and free proton fraction Yp for the accretion onto 
a non-rotating black hole, a = 0, with mass Mbh = 3Mq. The 
numbers of the label 1,2,3,4,5 denote the cases for accretion rates, M = 
O.O5,O.2,O.5,O.7,1(M0 sec"!). The radius r is normalized by the gravita- 
tional radius . 
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Figure 7. The same as Fig. 5, but for the cases of less accretion rates with 
a = 0.9 and Mbh = 3Mq. The numbers of the label 1,2,3,4,5 denote the 
cases for accretion rates, M = O.O()5,O.O1,().O2,O.O5,O.1(M0 sec"!). 
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Figure 8. The same as Fig. 6, but for a rotating black hole with a = 0.9. 
The numbers of the label 1,2,3,4,5 are same as in Fig. 7. 



Figure 6. The profiles of of electron neutrino Y^^ , electron anti-neutrino 
Yv^ , heavy neutrino Y^^ . Others are same as Figlsl 



?)7ip = {fJ.n — fJ,p)/kT, i.e., Un/up « exp(r;„p). The locus of 
accreting matter in the p — T plane is traced to the increase of 
chemical potential of fi„ — Hp (see Fig. 9). The larger accretion 
rate moves the locus into the larger potential in ij.„ — fip. In the typ- 
ical case of the accretion its factor rjnp increases from 5 at the outer 
boundary of the disk to 12 at the inner side of the disk (see Fig. 
10). The accretion flow with large efficiency of neutrino loss pro- 
duces the locus with relatively low temperature in the p — T plane, 
T « 2 ~ 5 X 10^'^ K, which means that the larger accretion flow 
has the larger number ratio, n„/np. These larger number ratios, 
rin/rip ~ 10'^^*, result in the larger differences in the composition 
and opacity between i/e and ■ 

The high dense matter and the sufficient time of reaction pro- 
ceed the electron capture on a proton. The electron fraction Ye grad- 
ually decreases along the flow. The fraction of electron neutrinos, 
Yl/^ , rapidly decreases along the accretion flow, while the fraction 
of antielectron neutrinos, Yi^^, gradually increases. These fraction 
changes are brought about mainly through the following reactions: 



e + p 



^ Ue +n, 
e ■ + n ^ i>e + p, 
e + e'^ ^ Ue + l^e 



(55) 
(56) 
(57) 



The fraction of heavy neutrinos Y^^ increases along the flow due to 



the reactions, e + Vi + i>i and N + N ^ N + N + Vi + i/i. 

In some cases, rh ^ 0.05 for a = 0.9, it exceeds that of Vc at the 
inner side of disk, Y^^ ^ Y„^ . The neutralization of the matter is 
remarkable, Ye ^ 0.03 for m ^ 0.02, near the inside boundary, 
r ~ Tms, where the equation of state for dense, hot matter drasti- 
cally changes. The composition of the accreting matter then rapidly 
changes. 



4.3 Emissivities and opacities of neutrinos 

Let's evaluate the emissivities of neutrinos produced by the ac- 
cretion. The emissivities by various types of reactions are shown 
in Fig. 11, in the typical case of the accretion. The most effi- 
cient emissivity is produced by the positron capture on neutron, 
+ n ^ i/e + p. At the outer side, r ^ 20rm3, the pair produc- 
tion, e + ^ i^e + Vc and the nuclei reaction, e + ^4 — > i/e + 
also contribute to the total emissivity. At the inner side of the 
disk the URCA process and nucleon — nucleon bremstrahlung be- 
come efficient. Though the emissivity by antielectron neutrinos 
qv^ is dominant, the other type of neutrinos, Vf^,Vx, also pro- 
duce some efficient emissivities. The emissivity by heavy neutri- 
nos, g^^, becomes comparable with q^^ near the inner boundary, 
qa^ ~ qv^ ~ 10^*(erg/cm''/sec). The emissivity increases sharply 
along the flow and the total emissivity is expressed approximately 
as g = q^^ + q^^ + q^^ ^ 10^'^r~'^ '^^(eig/cm^/sec). 

The opacities of absorption and scattering for each type of 
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Figure 9. The loci of accreting gas for the accretion rate M = 
0.01, 0.1, 1(Mq sec-i) on the p(g cm" 3) — r(K) plane. The contours 
of chemical potential fin — /J-p with constant electron fraction Ye = 0.3 
are also plotted. 
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Figure 10. The profiles of degeneracy factors for rjnp and rie . The profiles 
denoted by the accretion rates rh=0.1, 0.05, O.O2(M0sec~^) are shown, 
where a = 0.9, and those by m = l(MQsec~^) are in the case with 
a = 0. The mass of a black hole is set to be Mgu = 3Mq. 
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Figure 11. The profile of total emissivity of neutrinos q and those of anti- 
electron neutrinos q^^, electron neutrinos g^^, heavy neutrinos q^j. are 
plotted in the case of the accretion with Mbh = 3Mq, a = 0.9, and M = 
O-IMq/sbc. The emissivities by /3 process qp, positron capture on neu- 
tron e~ , e+ pair annihilation gee, URCA process quRCA^ nucleon- 
nucleon bremsstrahlung qnn, plasma process qpiasma^ and nucleus reac- 
tion qji are also depicted. The unit of emissivity q is (erg cm^'^sec^^). 
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Figure 12. The profiles of absorption opacity Kaii^i), scattering opacity 
Ksi^i), effective optical depth Teff(yi), and total optical depth r^^, for 
Vi = Ve,Ve,i'x, are plotted. The neutrino-nucleus opacity is labeled by 
KaA- The unit of opacity k is (cm~^). The parameters of accretion are 
same as in Fig. 1 1 . 



neutrinos are shown in Fig. 12. The absorption opacity of Ve. is very 
much larger than those of i^e, J^i^; i^a{ve) ^ i^ai^e) ^ f^ai^'x)- 
The electron neutrino is absorbed mainly by dense free neutrons 
through the reaction, Ve + n ^ e -f p, while the antielectron neu- 
trino Ve is little absorbed by rare free protons, Ve+P — » e++n. The 
scattering opacities of Ve and j/e are comparable except near the in- 
side boundary where the chemical potential of u^. is to some extent 
larger than that of Ue, rja^ > rj,^^ , which brings about the larger 
scattering opacity of i/e(see Table A2). The both total optical depths 
of and Ve, are very large with same order of magnitude at the in- 
side of the disk, Tj^^ — Ta^a + Ta^s > t^,^ — r^^a + r^^a ~ 10^. 
However, the difference of effective optical depth between r^j_^ 
and T^l^ is remarkable, i.e., r^f^C^ 30) > r^f^i^ 1), where 
the effective optical depth is defined by tI;^-^ = ^/tI^^t^.. Most 
of i/e are absorbed in the disk and their energy is deposited in the 
accreting matter. On the other hand the antielectron neutrinos i/e 
cany out the thermal energy from the disk to the outer atmosphere. 
The effective optical depth of i/e is approximately described as 
rl^y 6 X 10"^ T^^ {m)~^m. It becomes over unity at the radius, 
r 14r9(MB/^/3M0)-^/^(A^/(O.lM0/sec))^''^ where the 
electron neutrino sphere(disk) is formed. The sphere of antielectron 
neutrinos or that of heavy neutrinos are not formed when rh ^ 0.1. 
The total optical depths of Ve and become over unity, r^^ , Tj^^ ^ 
1, where r < 24rg{MBH /SMq)-'^''^{M /{O.IMq /sec)^^^ . 

The flux densities of neutrinos, F,^. (r), are shown in Fig. 11. 
The flux density of electron neutrinos, F^^, doesn't increase along 
the accreting flow. The ratio of the flux density of i^e to that of 
becomes large to be Fa^/Fu^ ~ lO^^** at the inside of the disk. 
The very large difference between F^^ and seriously restricts 
the efficiency of the annihiration of and at the outside of the 
disk, which results in the di fficulty in the formation of a fire ball by 
the annihiration of Ve,Ve jPopham. Wooslev & Frver' 1999). The 
flux density of i/^ rapidly increases near the inner boundary and 
reaches to be comparable with Fa^. The profile of the total flux 
density is expressed approximately as F{r) — F^^ + Fa^ + F^^ ~ 
3.8 X 10'"V^2-^(erg/cmVsec). 



4.4 Luminosity and mean energy of neutrinos 

The angular momentum of a black hole changes the efficiency of 
viscous heating through the shear stress ar^. The energy equation 
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Figure 13. The profiles of total flux density of neutrinos F, flux density of 
antielectron neutrinos Fi,^ , that of heavy neutrinos F^^ , and that of elec- 



tron neutrinos _F^^ . The unit of flux density F is (erg cm 
parameters of accretion are same as in Fig. 1 1 . 
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Figure 14. The profiles of relative 
rates of heating and cooling. The rel- 
ative heating rates due to viscosity 
fvis, due to compression fadv and 
the relative cooling rate due to neu- 
trino loss fi, are plotted in the case 
witii a = 0, Mgu = 3Mq , M = 
O.IMq sec-l. 



Figure IS. The same as Fig ll4l but 

for the case of a rapidly rotating 
black hole with a = 1, Mgu = 
3Mq are depicted, where M = 
O.O2M0 sec-l. 



is expressed such as 



dE 

at 



= O.S 

= 0, 



X 10''^(erg/sec/cm^)r 



(58) 



where f^is , fadv and are the factors of energy flux by vis- 
cous heating, by advectional heating and by neutrino cooling. 
These factors, fvis, fadv, fi^, for a — and 1 are shown in 
Fig. 14 and 15. The absolute values of f^i^ and f„ are orders 
of unity except for the inner side of the disk with a = 1. In 
the case of extreme Kerr, a = 1, the efficiency of viscous heat- 
ing fjjis decreases to be /„is « 0.05 at the inner boundary of 
the disk. The factor of advection, fadv, is relatively small and 
positive, fadv ~ 0.05. The accretion with the Keplerian angular 
momentum heats up the flow by the compression. On the other 
hand the advection dominant accreting flow(ADAF) cools the flow 
by the extension of the fl uid ( Di Matteo. Perna. & Naravan 2002 ; 
iNaravan. Piran. & Kumar; 2001 ; Pooham. Woos lev & Frver..l999.) . 

If the luminos ity of each flavor of neutrino and its mean en- 
ergy are observed ( Halzen & Hoooei 2002; McDonald A.B. et alJ 
[2003; McKinsev & Coaklev 2004), the physical structure around a 
black hole will be made clear. Though the electron neutrinos, j/g, 
are somewhat absorbed in the accreting matter at the outside of the 
disk, the antielectron neutrinos, i/e, and heavy neutrinos, Ux, are 
little absorbed there. The energy flux of i/e at the surface is very 
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Figure 16. The neutrino luminosity L(erg/sec) versus mean energy of 
ejected neutrinos i?(MeV) observed at infinity. The luminosity L and mean 
energy E for the variation of specific angular momentum of a black hole, 
a = — > 1, are depicted by the dot-dashed line labeled as m = 3[a = 
0, 1], where the other parameters of accretion are fixed to be Mgu = 
3M0,M = O.lM0/sec. Those for M_g// = 1, 2, 3, 4, 5, 6, 8, 1O(M0) 
are connected by the line labeled as a = 0.9[m = 10, 1], where a = 
0.9, M = O.lM0/sec. Those for M = 0.02, 0.05, O.1(M0 sec" i) are 
labeled as a = 0.9[dotm = 0.02, 0.1], where a = 0.9, Mbh = 3Mq. 
Those for a = 0, 0.7 are similarly shown. 



small in comparison with those of other types of neutrinos. Thus 
we calculate the luminosity L^. and the mean energy E,^. for each 
type of neutrino Vi emitted from the surface of the disk, and then 
discuss the effects of the absorption of at the outer atmosphere 
of the disk if necessary. The luminosity measured at infinity is ex- 
pressed in the Boyer-Lindquist coordinate as Li,. = J q^. {r)dSz, 
where gj. (r) is the z component of the flux density vector at the 
radius r, and dSz is the z component of the surface element vector. 
The component gj. is connected with the local surface density F^- 
such as gj. — e\Fy^, where e\ is the t component of the orbiting 
basis (15* <23> . Thus the luminosity is expressed by 



7(1 + v^V(^^))F^^rdr. 



(59) 



The mean energy E^^ is defined as follows. When the neu- 
trinos are thermally equilibrium in the matter, the local mean en- 
ergy of neutrinos e^^ (r) is expressed by the matter temperature 
T(r) and the degeneracy factor of neutrino rj^^ir), e^i{r) = 
kT{r)F3{r]^-{r)) / F2{r]„.{r)), where Fk{rii,^) is the Fermi inte- 
gral Fk{ri) = {1 + ex.p{x — T)))'^ dx. The number per 
unit time emitted from the differential ring over the width of ra- 
dius r ~ r + dr is dNv^{r) = dF^^{r) /i,^^{r), where dF^.{r) 
is the energy flux of neutrino from the ring. The total number of 
neutrino is Ni,- — dN„-{r). Thus the mean energy is defined 
hyE^. = L^JN^.. 

The total luminosity for all types of neutrinos is — 
y^^ . Lv^ . The mean energy in all types of neutrinos is simply de- 
fined by Ejj — Li, / ^2 ■ ^vi- We set the quantities of neutrinos 
within Tms to be same as the values at the boundary r = rms- The 
maximum radius of the integral is taken to be Tmax = lO^rms- 

The total luminosities and the mean energies measured at in- 
finity, Lv, El,, are plotted in Fig. 16 for some typical black holes 
with m = 1 ~ 10 and a = 1 and for the accretion rate m = 
0.02 ~ 1 . The rotating black hole has a small radius of the inner 
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boundary, r^s, where the high density of accreting matter increases 
the emissivity of neutrinos, which brings about the large luminosity 
with high mean energy Ev It is superior to the inefficiency of 
viscosity even if a ^ 1 and f^is 0.05. When the mass of a black 
hole is Mbh = SM© and the accretion rate is M = O.IM0 /sec, 
the ratio of luminosity is L^(a = 1)/Ly(a = 0) « 3.5, and that 
of mean energy is Ev{a = 1) /Ev{a = 0) ~ 2.5. The luminosity 
is proportional to the accretion rate, oc M, but is independent 
of the mass scale of a black hole, Mbh, that is, the constant ratio 
of the gravitational binding energy of accreting matter is released 
into the neutrino's energy through the viscous heating. Its ratio in- 
creases according to a. The mean energy Ei, is proportional to M 
but is inversely proportional to Mbh, Ev oc M« »M5^"^ When 
each type of neutrino ejected from a GRB is observed, the physics 
around a black hole will be made more clear. The observed lumi- 
nosities of neutrinos and their mean energies will determine the 
physical parameters of a central black hole, a, Mbh, and the ac- 
cretion rate, M, from the diagram shown in Fig. 16. 
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Figure 17. The variation of net heating rate SQ due to the increase of 
temperature ST within the diffusion time. The variation of the rate by the 
viscous and compressing heating <5Q+, and that by the neutrino cooling 
are plotted by dashed lines. The scattering optical depth of antielec- 
tron neutrinos Ts is plotted by dot-dashed line. The parameters of accretion 
are M = O.SMQsec"^, Mbh = SM©, a = 0.9, = 0.1. 



5 DISCUSSION 



We have studied the accretion in simplified model in which the disk 
is assumed to be stationary and its vertical structure is treated by 
one-zone model with the thickness, H. These simplifications might 
restrict the properties of the accretion disk. We first discuss the ther- 
mal stability of the disk. In the stationary state, the cooling rate due 
to neutrino's loss must be equal to the production rate of thermal en- 
ergy in the disk even if neutrinos are scattered in many times within 
the disk. The accretion disk is thermally stable since the antielec- 
tron neutrinos i^e are little absorbed in the disk and can sufficiently 
carry out the perturbed thermal energy. However in the short dura- 
tion within the diffusion time, tdiff ~ TsH/c ~ lO(ms), the per- 
turbed energy is stayed in the disk. The heating or cooling time then 
becomes theat = tcooi =~ e/i? ~a few ms, i.e., theat < tdiff- 
If the accretion disk is thermally unstable, the perturbations rapidly 
grow. This stability should be decided by the fully nonstationary 
treatment of accretion disk. Here we discuss it in brief by the pertur- 
bation analysis. When the scattering optical depth is large Ts 2> 1 
and the absorbing optical depth is negligible, Ta ^ 1, the neutri- 
nos stay in the disk in times of the crossing time t cross = H /c. 
When Ts <C 1, the perturbation of the cooling rate is expressed as 
5Q^ — SFv — 5{q^H). We provide the perturbation of the cool- 
ing rate in the thick case as SQ^ — 5{q^H/Ts). The perturbed 
cooling rates by i/e and i/x are then set to be SQ^^ + SQ~^ = 
5{qi,^H/{l + Ts(ve))) + S{q^^H/{l + Ts(ux))). Since the elec- 
tric neutrinos u^. are almost absorbed, the perturbed cooling rate is 
expressed as 5Q~^ = 5(g^^i//(1.5r„^ar„^ + \/Zri,^a + 1)). 

The perturbation of the heating rate is brought about by the 
viscous heating and the compression, SQ^ — SQvis + SQadv We 
thus evaluate the thermal stability for the perturbation of tempera- 
ture ST with the restriction of constant surface density E, 



91ogQ+ 
91ogT 



+ 



dlosT 



(60) 



The results for SQ are shown in Fig. 17, where m = 0.5, m — 
3, a = O.Qandayis = 0.1. The sign of SQ becomes positive when 
Ts{vc) ^ 1. The sign of 5Q^ = (c'log log T)e also be- 
comes positive when rs(i'e) ^ 10, since the opacity of h'e increases 
according to the rise in temperature, ST > 0. For the region where 
the scattering optical depth of Ve is larger than unity, rs(i/e) ^ 1, 



the accretion disk becomes thermally unstable at least within the 
duration of the diffusion time. 

We next discuss the interaction between neutrinos and nu- 
cleons. The spectrums of neutrino are assumed here to be Ferm- 
Dirac distributions with local temperature being equal to that of 
nucleon medium. In the neutrino sphere(disk) the neutrinos are 
thermally equilibrium with nucleons but in the outer layer above 
it the above assumption is in general not consistent. In the outer 
layer the spectral temperature defined by the mean energy of neu- 
trinos, T, = EvI'i, becomes higher than the nucleon's tempera- 
ture iRaffelt 2001). Since the mean energy of neutrino exceeds the 
rest mass energy of electron, E^ > ra^t? , the collision between 
neutrinos and electrons effectively transfer the momentum and en- 
ergy from neutrinos to electrons I Thompson, B urrows. & Horvatl] 
^000). The dynamic structure functions which include the full kine- 
matics of -nucleon scattering produce also the tranfers of energy 
and momentum between neutrinos Vi and nucleons. In the diffu- 
sion zone in which the neutrinos many scatter the electrons and 
nucleons, Ts ^ 1, the convection originated by neutrinos may 
work important effects on the dynamics of the medium. In our 
accretion model the effective optical depth, t^] f = yfr^s, ex- 
ceeds unity only for v^- The thermal equilibrium disk of is 
formed at r ~ lOrg. The scattering optical depths of all flavors, 
Ts(ye),Ts(pe),Ta(yx\ excced unity. In the diffusion zone above 
the disk, z ^ H, the neutrinos might push the nucleon matter out- 
ward and might drive the convection or wind or jets. 

Our results for the hot, dense matter accretion onto a black 
hole indicate that the emissivity of electron antineutrinos q^^ is 
very efficient and the fraction of free protons is extremely minute, 
which allows the accretion disk effectively cooled by u^.- On the 
other hand the electron neutrinos are almost trapped by many 
free neutrons. The accretion disk cooled by neutrinos could be ther- 
mally unstable at least within the diffusion time, tdiff ~ lO(ms). 
Our treatment of the accretion is very simple, stationary model with 
the homogeneous disk in the vertical direction. The compactness 
problem for GRBs requires the formation of relativistic jets. The 
accretion should be further solved in unstationary, 2- or 3- dimen- 
sional treatment with the transfers of neutrinos. 
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APPENDIX A: REACTION RATES OF NEUTRINOS 

The neutrino cross sections, opacities, emissivities and reaction 
rates used in our model are listed in Table A 1-4. The numerical 
functions used in the reaction rates are listed in Table A5. 



APPENDIX B: TRANSFERS OF NEUTRINOS IN A 
PLANE-PARALLEL DISK 

The analytic formulas of neutrino's transfers are derived by 
using the simplified treatment of radiative transfer given by 
IPopham. & Naravaril <T995I) . Introducing the two streams, with in- 
tensities Z"*", I~ moving at an angle cos 6 — l/y/S relative to +z 
and —z, define the net flux density j and the to tal intensity j such as 
j = (7+-/-)/2V3and j = (1+ + J") /2 I'Popha m. & NaravarJ 
1995"). The transfers of neutrinos are then described by the produc- 
tion rate of neutrinos, h'^ , and the absorption opacity, scattering 
opacity and the total opacity, Ha, and k — Ka + k.^, as 



di 

dz 

§_ 

dz 



-Kaj + h 
-3K.j. 



(Bl) 
(B2) 



The boundary condition at the center of the disk is j = 0. To 
derive the simplest solution, let us assume that (n+ — Kaj), k are 
independent of z tiHubenvtl990i') . Equations <BH and <B2> can then 
be integrated to give 



j{z) = (n^ - Kajc)z 

= jc - T;li{nt - Kajc)z'^, 



(B3) 
(B4) 



where n J and jc are the values of h'^ and j at z — 0. The boundary 
condition at the surface of the disk is j~ = 0, i.e., j — \/3j- 
Applying this boundary condition we find 



V3 + 1.5T . + ^ 
Jc = ]= Uc a, 

1 + VSra + 1.5rTa 



(B5) 



where Ta is the absorbing optical depth, Ta = KaH, and r is the 
total optical depth, r — kH. The escaping flux density of any sort 
of neutrino, fi, is then given by 



= j{H) = 



ntH 



(B6) 



where r,^; and Ti,^a are the total optical depth and the absorbing one 
of Ui, and n^. is the production rate of Vi per unit volume. 
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Table Al. Cross section 
No. Neutrino Reaction 


Cross Section 


1 


+ e + p 




2 


i'e + P — > e+ + n 




3 


Vi + p Vi + p 






Vi + n Vi + n 




5 


Vi + e-* Ui + e 


= <To| (Ci 



{Cv+Ca)^ + UCv-Ca? 



°' Tliese cross sections are given by Burrows(2001). A convenient reference neutrino cross section is cfq 



*cin'^. gA is the axial- vector coupling constant (~ —1.26), 9\y is tlie Weinberg angle and sin^ 



7r(hc/27r)4 



~ 1.705 X 



10- 

1.29332MeV, and for a collision in which the electron gets all of the kinetic energy 

correction and is approximately equal to (1 + l.l£i/g/mnC^). Cv = 1/2 + 2 sin^ for ^^e and Pe neutrinos, Cv 
for Vx and Vx neutrinos, Ca = +1/2 for t-e, S'fi and Vr neutrinos, and Ca = —1/2 for Pe, ^"/^ and Vt neutrinos. 



0.23, Q = ninc^ — nipc^ 



Q. Wm is the weak magnetism/recoil 
-l/2 + 2sin2 Ow 



Table A2. Opacity 

No. Neutrino Reaction 



Opacity 



References 



Ue + n ^ e + p 



Ue + p - 



- N + N - 



v, + N - 
i>i + A- 



Vi+e 
*Vi + A 



N + N 



J d''p/..e(e)<T;^,n(e)[l-/e(e 



: 4.3 X 10-6(cm-l)pi2y„pT2j^S^J^ 



+Q)] 



/3 
a 

NN 



4.3 X 10-6(crn-i)pi2ypnT2j 



F,F.. 



1.43 X 10-9(cm-i)(X„pi2)^T,2i^{^( 2+o.862T„ ) 



1/2 



^ ^aoUAei^-^p-" jNp{Z)Nh{N) 
1.2 X 10-6(cm-l)pi2yAT2j/A(r;) 



Kj, ^ 1.1 X 10-6(cm-l)pi2T2^C,,Aryjviv^{^ 

Kl « 1.66 X 10-6(cm-i)yePi2T2^^{g 



-faaiiAA'^ 



kT 

me c. 



■£4(37^ 



; 2.3 X 10-'^(cm-i)pi2r2^A2y^ 



Burrows(2001) 

Burrows(2001) 
Raftelt(2001) 

Bruenn(1985) 

Burrows(2001) 
Burrows(2001) 

Bruenn(1985) 



The absorption opacity Ka and scattering opacity Ks are expressed with unit cm ^ . Fk (rf) is the Fermi integral M = f "" ^ <^^^ 

J (J i-|-cxp l^a; Tj J 

The complex integrals, F^{ri^^ , rie,q) and ~ (17^^ , ?7e), are shown in Table A5. The doubly indexed quantity Y„p is the number fraction 
defined as Ynp = Vnp/nB, and Yp„ is Yp„ = exp (r/p — r]„)Ynp, where rjnp is the dynamic structure factor for neutrino-neucleon in- 
teraction, = / ^/„(e)[l - fp{e)] = expr>7p-"^")-i - ^^"^ '^^^ = l+|ml"(„.^,,(l) - = 1 if = n 
and C's,N = 0.827 if N = p. A is the mean mass number of heavy nuclei, ua is the number density of heavy nuclei and Ya is the 

fraction, Ya = UA/ns- is the mass fraction of neutrons. fA{v) is fAiv) = e^"~'^p~'' jNp(Z)Nh(N) p^^^° '^"^ where 
q' PS 7j„ - 77p + A, A = 2^^. Np{Z) and Ar^(Ar) are Arp(Z) = {0 for Z < 20; Z - 20 for 20 < Z < 28; 8 for Z > 28}, and 
Nfi{N) = {6 for AT < 34; 40 - AT for 34 < AT < 40; for N > 40}. The density pi2 and temperature Tn are normalized vales by 
10i2(g/cm3)andlO"K 

APPENDIX C: CHANGING RATES OF POSITRONS „ _ 1.8/ kT f/„ w ^ i f(„ w 

ALONG THE FLOW " ^ \h^) < ' " ^^'^^^ ' " •^^''^^ ^ > ' ^^'^ 



The number density of positrons which are equilibrium with ther- (Landau and Lifshitz 1964), where the brackets indicate the mean 
mal photons is described as values of the phase space blocking factors. These factors are eval- 
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Table A3. Emissivity 
No. Neutrino Reaction 



Emission Rate 



References 



N - 



hn ■ 

■ N - 



■ n - 



■V 
N - 



■ n - 



■ Vi+Vi 



Si 7.4 X 1033(ergcm-3sec-i)pi2T-f^yp„F5"(r?^,,r)e,g) 

q-^{v^) ^ 7.4 X 1033(ergcm-3sec-l)pi2rfir„pFg- (r;^„,»?^+,g) 

qNN{vi,Vi) ^ 7.3 X 1030(ergcm-3sec-l)(XArpi2)2T5^5 



Si 2.06 X 1033(ergcm-3sec-i)pi2yAf Afp(^)AfhWTfiF5' ~{-n^,,Ve) 

quRCA(ye) = 2.7 X 1034(ergcm-3sec-i)pif TfiFpCrye.rj^J 
q~,{ue,Pe) ~ 1.20 X 1033(ergcm-3sec-i)T9FQ^(»7i,J 
q-y{u^) X 8.3 X 1030(ergcm-3sec-l)T9FQ.^(»?,^^) 
qeeiiye,Pe) ^ 1-24 X 10^3 (ergcm-Ssec" ^ )T9^ 
qe4ux,Px) « 5.49 X 1032((,rgcm-3sec-i)Ti\/^J»?) 



Bun-ows(2001) 

Burrows(200I) 
Bun-ows(2001) 

Bruenn(1985) 



Shapiro et al.(1983) 
Ruffertetal.(1996) 
Ruffertetal.(1996) 
Thompson et al.(2000) 
Thompson et al.(2000) 



The emission rates q are expressed with unit erg/cm'^ /sec. The complex integrals, F~ {rju^,rie,q), F~^{rii)^,ri^+ ,q) and P'j? {rju^), 
are shown in Table A5. Fq and f^^ for i = e,x are FQ('qe,'ri,y^) =< 1 — fe{se) >< 1 — fvei^t^c) > fi^iiv) = 

2F4(0)F3(0) < >ee< l-/>^i >ee-i^Q7 IS FQT,(?yi, J =7'=e ^^({I + Ct+I) ) < l-^eC^) >7< 

1 — fv^{e) >7, where 7 and 70 are 7 sd 70 -^/""^ /3 + r)| , 70 = hQ,Q / = 5.565 X 10~^. The mean blocking factors 
< 1 — fu^ (e) >ee, •• are described in Table A5. The density pi2 and temperature Tn are normalized vales by 10^^(g/cm'') and lO^^K 



Table A4. Fraction Change 
No. Neutrino Reaction 



Fraction Change 



1 



+ n ^ e + p 



A^e + A' 



n + n—>n + p + e-\-Ue 



N + N - 
e + e+ ^ 



Y„pU,{e)[l -fe{e + Q)] - Yp„Me + Q)[l - UAe) 
Si 9.0 X 102(sec-i)Tii [Y„pF+{7j^^ ,Vc,q)- Yp„F-{r,^^ , r)^,q)] 

2.4 X 1 f^*^ — V ,2lAT I 7\AT, ( M\ I aVn—V-D-Q' T?1 ' 



+ n 

^ N + N + Vi + Ui Y^,{NN) 

- Vi + Vi 



102 {&ec-^)Tl^YAljNp{Z)Nn (N) (^e"' 
Y4URCA) Si 

Si 3.3 X 103{scc-1)p,-'/'t,7i^{^ < 1 - Ms) >„„< 1 - > 
i;+(Pep) Si 9.0 X 102(sec-i)Tii[yp„F+(j7s,,»7^+, -g) - Y„pF^ [rji,^ ,r]^+ , -q)] 

2e„.ni5 ~ J--Jl>'t-'- jvs.^pi2-iij^ 1 2+0.862Tii 

y.Jee+) Si 36.3(sec-i)Tf,/3,,J»)) 
y.Jee+) Si 31.2(sec-l)Tf,/3,., (r,) 
y.,{7) « 1-4 X 102(sec-l)F7(r,)p-iTf, 
y.j7) ~ 1.0{sec-i)F7(r,)p-irf, 



1/2 



The parameter 6 is b = 47r J . The complex integrals, F^ (Vt^etVe) and ^4 (^i^o^e), are shown in Table A5. /3,iyj(r;) and 
FQ.,(r,,J are/3,,^(r,) = (F3(0))-2 |F3(77,)F3(r?^+ ) < 1 - M{s) 1 ~ -i^3(r?.. )i^3('?s. ) < 1 - fe{e) 

1 - /e+(e) >'^ii^, |. and FQ^,{r)^J = ■y^e-'y^{{l + (7 + 1)^) < 1 - /^,(e) >7< 1 - /i;,(e) >7. The mean blocking factors 



. are described in Table 5. Xn is the mass fraction of neucleons. 



uated with the average energy of electrons and positro ns produced 
by e~e+-pair creation fRuffert. Janka. & SchafeJl99d) : 

< 1 -. f ) >={l+cxp[-(e+- 7,,+ )]} \ (C2) 

The average energies, e"*" and e~ , are approximately given by 



77r^ 
180C(3) 



kT, 



(C3) 



where 6^+ and 11^+ are the energy and number densities of 
positrons. When the fraction of positrons becomes large, Ye+ ^ 
0.01, the changing rate of positron number, h^+, due to the reac- 
tions, e'^ + n ^ Us+p , e + e+ ^ i^i + i/i, cannot be negligible. 
We calculate the changing rate of the positrons along the flow by 
using the following equation: 



-1 r-r-. 



(C4) 
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Table AS. Integrals and Blocking Factors 
No. Integrals and Blocking Factors 



-''fc Wi',»?e,g) - Jq i + cxp (.J- - /eig + aijjda; - (i+oxp (a;-,,„))(l+exp (,,,-9-2,)) 

(£)>nn«{l+Cxp[^.-^(^]} ' 



- fe 
< 1-/. 



E + — ^7e)){l+exp {rjy — a:)) 
{x + q' — 77e) + l)(exp (77i^ — a::) + l) 
l + oxp (a;-77„)){l+exp (ric—q'—x)) 



< 1 - /,^(e) >„„ai |l + cxp [77^ -pT- 

< 1 - Uiis) >ee~ |l +exp 



£5(^1 \ ^ 



1 ( F5('?e) 



< 1 - fe{e + g') >~ |l + cxp [ije 



f5('7.e) 



< 1 



/-^.(e) >7~ (l + exp [^,, -(l + i^ll)]) 



The quantity 0(x) is the unit step function, 0{x) = if x < 0, and ©(a;) = 1 if x ^ 0. 



where is the net production rate of positrons, — —hjj — 



APPENDIX D: DYNAMICAL STRUCTURE OF THE DISK 

Dynamically equilibrium structure of the disk is given by the Euler 
equation. The equation of motion of the fluid, h^T^,^ — 0, for the 
steady, axially symmetric and rotating ideal fluid can be expressed 
in the total differential form: 

— — — = ^'^{—dv + vf^-^dip — v^^-)V^dlnuj) = dU. (Dl) 



The equipressure surfaces are given by the equation U ^constant. 
The quantity U — U{p) is equal in the Newtonian limit to the total 
potent ial (gravitational plus centrifugal ones) e xpressed in the units 
of iAbramowicz. Jaroszviiski Sikoralll978h . We solve the disk 
structure in terms of cylindrical coordinates ro, 2, and (p instead of 
radial coordinates which are combined through the relations, 
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If no pressure force exerts in the to direction, the equation ijDlfl 
determines the velocity of the geodesic circular orbit, v^^y Let's 
introduce the new variables, a, m, s, c defined by a = a/r, rh = 
M/r, s = sin 9, c = cos 9. The following functions are defined 
by: 
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The circular velocity is then given by 



-B^ ± VB'i + A^C^ 
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This velocity at the equatorial plane is expresses as 
ei(r) „ _ 1 
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The vertical force is described by the gradient of the equipo- 
tential, Vzf/. The first order of the expansion of dU/dz in z/r{<^ 
1) gives the approximate expression such as 
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